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ABSTRACT 

We calculate the electron acceleration in random superluminal strong waves 
(SLSWs) and radiation from them by using numerical methods in the context 
of the termination shock of the pulsar wind nebulae. We pursue the electrons 
by solving the equation of motion in the analytically expressed electromagnetic 
turbulences. These consist of primary SLSW and isotropically distributed sec¬ 
ondary electromagnetic waves. Under the dominance of the secondary waves, all 
electrons gain nearly equal energy. On the other hand, when the primary wave 
is dominant, selective acceleration occurs. The phase of the primary wave felt by 
the electrons moving nearly along the wavevector changes very slowly compared 
to the oscillation of the wave, which is called ’’phase locked”, and such electrons 
are continuously accelerated. This acceleration by SLSWs may play a crucial role 
in the pre-acceleration for the shock acceleration. In general, the radiation from 
the phase-locked population is different from the synchro-Compton radiation. 
However, when the amplitude of the secondary waves is not extremely weaker 
than that of the primary wave, the typical frequency can be estimated from the 
synchro-Compton theory by using the secondary waves. The primary wave does 
not contribute to the radiation, because the SLSW accelerates electrons almost 
linearly. This radiation can be observed as a radio knot at the upstream of the 
termination shock of the pulsar wind nebulae without counter parts in higher 
frequency range. 

Subject headings: pulsar wind nebula: general — radiation mechanisms: general 
— particle acceleration 


1. Introduction 

Particle acceleration is one of the most important physical processes in the astrophysics 
and plasma physics. The acceleration by the strong electromagnetic waves (EM waves) has 


been investigated for a long period in the field of laser physics (e.g. Jory and Trivelpiece 
1968, Sarachik and Schappert 1970, Karimdabadi et ah 1990, Kuznetsov 2014). In the held 
of astrophysics, such mechanisms were investigated mainly in the context of the pulsar. The 
magnetic dipole radiation was thought to be emitted from the pulsar which has inclined 
magnetic axis to the rotation axis (Pacini 1968). This EM wave was regarded as a ’’strong 
wave” in the vicinity of the pulsar. The strength of the EM wave is generally dehned by the 
strength parameter a = eE/mcu, where E and u are the amplitude and frequency of the EM 
wave, m and e are the mass and charge of the electron, and c is the speed of light in vacuum. 
It is estimated to be much larger than unity for the EM wave around the pulsar, which means 
that this EM wave is capable to accelerate the electrons to a relativistic energy. Electrons 
dropped in the strong EM wave at rest are strongly accelerated toward the direction of the 
wavevector. The phase of the wave felt by the electrons changes very slowly (phase locking), 
since electron speed becomes close to speed of light. The strong wave is nearly a stationary 
EM held for these electrons, and they are continuously accelerated. If we assume the inhnite 
plane wave, the maximum Lorentz factor is 7max ~ not a (Gunn & Ostriker 1971). To 
achieve 7max ~ phase locking effect plays a crucial role. 

As is well known, the pulsar magnetosphere is not a vacuum as was assumed by the above 
papers, but is hlled with dense plasma (Goldreich & Julian 1969). Kegel (1971) pointed out 
that when the refraction index nr is signihcantly smaller than unity because of the existence 
of the plasma, the phase locking becomes inefficient since the phase velocity c/nr of the 
EM wave becomes signihcantly larger than c. Such waves are called superluminal strong 
waves (SLSWs). Moreover, if the plasma is sufficiently dense, even the propagation near 
the light cylinder is prohibited. When we consider the pair creation in the magnetosphere, 
the plasma around the light cylinder is much larger than the Goldreich-Julian density. The 
low frequency EM waves such as the dipole radiation from the pulsar cannot propagate in 
such a high density plasma. The magnetic energy is not carried by the EM wave but by the 
entropy wave in the magnetised plasma wind (i.e. striped wind in Goroniti 1990). 

Recently, the strong EM waves in the pulsar environment again attract many interests. 
It is pointed out that the entropy mode can be converted to the strong EM wave in the outer 
region of the pulsar wind (Arka & Kirk 2012). Moreover, it is numerically shown by using 
relativistic two-fluid simulation that this conversion can occur by the interaction with the 
termination shock (Amano & Kirk 2013). Such EM waves have superluminal phase velocity 
and their strength parameter a may be larger than 1. Around the termination shock of 
the pulsar wind, the electrons are accelerated by shock crossing and radiate synchrotron 
photons. If such SLSWs exist around the shock, they may affect the particle acceleration 
and radiation. 
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There are many unsolved problems for the particle acceleration around the termination 
shock of the pulsar wind nebulae. One is the injection problem. In general, to be injected 
into the shock crossing cycle, the particles have to be supra-thermal when they encounter 
the shock front. However, in the paper of Kennel & Coroniti (1984a), upstream plasma is 
assumed to be cold and all particles are accelerated. In other words, they assumed a very 
high injection rate. The lowest Lorentz factor 7 ~ 10® at the immediate downstream of the 
shock corresponds to the bulk Lorentz factor of the upstream wind T ~ 10®. These particles 
emit optical photons by the synchrotron radiation. Their model can explain the observed 
spectrum at frequency region higher than the optical range, but radio components were not 
discussed (Kennel Coroniti 1984b). There are some models for Crab Nebula which explain 
the radio components. For example, we have time-dependent one zone models (Tanaka & 
Takahara 2010, Bucciantini, Arons, & Amato 2011) and axisymmetric two-dimensional MHD 
(MagnetoHydroDynamics) models (Olmi et al. 2014, 2015). They reproduces the spectrum 
of the Crab Nebula by assuming the energy distribution of the electrons as a broken power 
law, but they do not specify the origin of the break Lorentz factor 7 ~ 10®. It may imply 
that there exists a pre-acceleration mechanism in the upstream of the termination shock. 
That is to say, we can consider a model in which the pre-acceleration mechanism makes the 
energy distribution broader, and only higher energy components are injected to the shock 
crossing cycle. If the SLSWs can make a non-thermal energy distribution in the upstream, 
it can be a good candidate of the pre-accelerator. 

Radiation from the electron in a strong (a > 1) EM wave is usually called synchro- 
Compton radiation (Rees 1971) or nonlinear inverse Compton (NIC) radiation (Gunn & 
Ostriker 1971). The deflection angle in one cycle of the electron motion in the interaction 
with a strong wave tends to be larger than I/ 7 . Asa result, the radiation signature resembles 
the synchrotron radiation. However, we do not know the radiation spectra when there are 
many waves. Moreover, even for the interaction with one SLSW, there are ambiguous points 
in the consideration of the typical frequency. To estimate the typical frequency of the 
radiation, we have to know the photon formation time (Akhiezer & Shul’ga 1987, Reville 
& Kirk 2010), which is the inverse of the cyclotron frequency mc/eB in the context of the 
synchrotron radiation. Previous studies considered the cases for which the phase locking 
effect is weak, which means that the velocity is oblique to the wavevector direction. For 
this case, the photon formation time is ~ mc/eE. In general, the motion of the electron in 
the strong wave can result in a different photon formation time. Therefore, the resultant 
radiation spectra can be different from the NIC or synchro-Compton theories. 

In this paper, we study the electron acceleration in the SLSWs. Also, the radiation 
signature in such situation is studied. We use numerical methods to investigate such highly 
nonlinear motions. The contents of this paper are as follows. In section [2l we discuss 
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the physical parameters for the computational study and describe the methods we use. The 
results are shown in section [3l Section IHpresents discussions including observational features. 
We hnally summarize this paper in section O 


2. Formulation 
2.1. Parameters 


Here, we estimate the physical parameters around the termination shock of the pulsar 
wind nebula by using the Crab Pulsar & its nebula as a representative system. First, 
we estimate the strength parameter a of the entropy mode in the immediate upstream 
of the termination shock by following the estimation by Kirk & Mochol (2011), in which 
they estimated a in the (Active Galactic nuclei) AGN jets. We note that this strength 
parameter is not identical to the strength parameter of the EM wave. The strength parameter 
of the entropy mode is defined by the wavelength Asw and magnetic field strength B as 
ei?Asw/27rmc^ in the observer frame. There are two feasible assumptions. One is that the 
magnetic field strength in pulsar wind is inversely proportional to the distance from the 
pulsar r as B (X r~^, (here we implicitly assumed the magnetic field configuration to be 
pure toroidal). The other is that the spindown luminosity of the pulsar Lgd is carried by the 
Poynting flux (high sigma) and to be isotropic. On the above assumptions, the entropy wave 
in the observer frame resembles EM wave due to the fact that electric field is perpendicular 
to the magnetic held, the strength ratio is nearly unity, and the wavevector is perpendicular 
to the electric held and magnetic held. We obtain a from the observed quantities as 


a = 



e^L 


sd 


m?c’ 


1/2 

~ 3.4 X 10^° 



/ Lsd 
^IQSSerg/s 


1/2 


( 1 ) 


where tlc = 1.6 x 10®cm is the light cylinder radius of the Crab Pulsar. The radius of the 
termination shock of Crab Nebula is ~ lO^rLc and spindown luminosity of the Crab Pulsar 
is ~ 6 X 10^®erg/s. Thus, the strength parameter of the entropy mode is estimated as a ~ 80 
at the termination shock. Since the EM wave is expected to be converted from this entropy 
wave, we can expect this wave to be SuperLuminal ’’Strong” Wave (SLSW). 

Next, we compare the inertial length and the wavelength of the entropy mode in up¬ 
stream and downstream by assuming that the typical length scale does not change even if 
some fraction of energy is converted from the EM energy to the kinetic energy, such as the 
magnetic reconnection (Lyubarsky & Kirk 2001). If the wavelength of the entropy mode is 
shorter than the inertial length, the MHD approximation breaks down and entropy wave can 




( 2 ) 


be converted to the other waves. The ratio of these length scales at the upstream is 

^ ~ 2.7 X 10^ 

^pjUp 

where T is the bulk Lorentz factor of the pulsar wind, = 10®cm is the wavelength of the 
striped wind in the observer frame, cUp^up = \/47r?7,upe^/m is the upstream plasma frequency, 
and riup is the comoving number density. Here we assumed that the plasma does not have 
relativistic temperature. On the other hand, constraint on Uup and T at the termination 
shock is obtained from the spindown luminosity. By expressing the spindown luminosity as 
a sum of isotropic kinetic and Poyniting fluxes, and substituting the radius of the termination 
shock of the Crab Nebula into it, we obtain 

r2(l + (T)nup = 2 X 10-2. (3) 

Here a is the ratio of the Poynting flux to the kinetic flux. Using this constraint, equation 
([2]) gives 

~ 3.8 X 10^ X (1 + (4) 

We can estimate at the immediate upstream of the termination shock by using a. At 
the light cylinder, the pair cascade models predict a 10^ (e.g. Hirotani 2006). It should 
be reduced by the magnetic reconnection in the wind region. If this process is extremely 
inefficient and a is close to 10^ at the termination shock, r]^p < 1 is realized. In this case, 
the entropy wave may break down by non-MHD effect. On the other hand, if we adopt more 
a conventional value of a < 10^, r]^p should be larger than unity. In this case, the entropy 
wave can survive in upstream. Hereafter we consider the latter case. 

The situation drastically changes in the downstream. First we estimate the ratio in 
the downstream r/down under the assumption that the entropy mode does not convert to the 
SLSWs. From this estimation, we can realize that this assumption is not appropriate. Then, 
we consider an alternative scenario. 



According to the observational fact, the bulk Lorentz factor of the downstream plasma 
is F ~ 1. The wavelength of the entropy wave changes due to shock compression, but their 
compression ratio measured in the downstream frame is at most 0(1). Therefore, this scale 
is ~ Asw in the downstream frame. On the other hand, the inertial length in the downstream 
frame increases drastically because of the thermalization. The ratio of them becomes 


A« 


hdown 


27rc/ CUp^down 




hup) 


(5) 


where 7th is the Lorentz factor of the thermal motion of the downstream plasma, Wp^down = 
^47r?7,downe2/7th??^ is the plasma frequency which takes into account the relativistic correc¬ 
tion, and ridown is the downstream number density. We assume 7 th ~ F, since the downstream 
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bulk velocity is nonrelativistic. As a result, ?7down is smaller than ? 7 up by a factor of F. This 
bulk Lorentz factor is estimated as 10^ < F < 10® by pair cascade models (e.g. Hirotani 
2006), many MHD models of Pulsar Wind Nebula (PWN), and the induced Compton scat¬ 
tering constraint (Tanaka & Takahara 2013). Using the constraint, we estimate ?7down to 
be 

10-® < r/down < 1. (6) 

Thus, MHD approximation is not adequate to describe the entropy waves in the downstream, 
and it should be dissipated (Petri & Lyubarsky 2007, Sironi & Spitkovsky 2011) or converted 
to some other waves (Amano & Kirk 2013). 

Before proceeding further, let us summarize the above estimates. The strength param¬ 
eter a of the SLSWs is estimated as 0(10). The wavelength of the entropy mode is expected 
to be longer than the inertial length in upstream, while the opposite is expected in the 
downstream. The expected values of their ratios are r|^^p > 1 and r^down ~ 10“^. 

The entropy mode cannot be sustained by the downstream plasma. It should be con¬ 
verted to EM waves by the interaction with the shock front and the EM waves propagate 
back to the upstream region (Amano & Kirk 2013). As we will see below, the entropy mode 
can convert to SLSWs even in the upstream. Hereafter we consider the propagation condi¬ 
tion of a SLSW, which is affected by the strength of the EM wave and the thermalization of 
the background plasma. 


First, we review the effect of the strength of the wave. The frequency of the SLSW 
should be ~ 27rc/Asw in the shock rest frame. Here, we assume that this SLSW propagates 
toward the shock front in this frame. The frequency of this SLSW in the upstream frame is 
27rc/PAsw For an ordinary (a < 1) EM wave in the plasma, the condition for the propagation 
is > Up. This condition coincides with ?7up < 1. However, for SLSW, the condition is 
modihed by the strength parameter. The dispersion relation of the SLSW for pair plasma is 
described by Kaw & Dawson (1970) as 


>2 

2 I ; 2 2 

u = , + k c , 


VT 




where k is the wavevector. Hence, the condition for the propagation becomes 

-2 2 


u 


> 


u. 


p vT+^ 


Therefore, the wave can propagate when a is sufhciently large even if ? 7 up > 1. 


(7) 

( 8 ) 


Next, we discuss the thermalization. The SLSWs in the over dense plasma (i.e. f] > 1) 
are unstable (Max & Parkins 1971), and also it is shown by the two-fluid simulation that these 
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waves generate EM waves and sound waves by the stimulated Brillouin scattering (Amano 
& Kirk 2013). As a result, the SLSWs thermalize the plasma in the upstream region by 
the dissipation of sound waves. If this mechanism works well, the upstream plasma gain a 
relativistic thermal energy. The plasma frequency becomes lower as Lorentz factor of the 
thermal motion becomes large as Up oc ' . Thanks to this nonlinear effect, the local 
plasma frequency in upstream becomes lower than the SLSW frequency. The instabilities of 
the SLSWs are not fully understood. The growth rate is only known for limited conditions 
(c.f. Max & Parkins 1972, Asseo et ah 1978, Lee & Lerche 1978). It tends to be low when 
a is not much larger than 1 and u is sufficiently larger than Up (Amano 2014). This is 
natural because the wave becomes an ordinary EM wave in vacuum in the limit of a 1 
and u 3> Up. 

As we have seen above, SLSWs can exist in the upstream region of the termination 
shock. In this paper, we study the electron acceleration in the upstream rest frame. We 
assume that a = 0(10) and u > Up. These assumptions are acceptable near the immediate 
upstream region. 


2.2. Setup 

We describe the EM turbulences by the superposition of the EM waves. In this paper, we 
ignore wave-wave interactions for simplicity. Such interaction may not be negligible and will 
be treated in future works. This description of the turbulence is based on Giacalone & Jokipii 
(1999), in which they calculated the transport of the cosmic rays in the static magnetic field. 
Recently, they studied particle accelerations using this scheme (e.g., Giacalone & Jokipii 
2009, Guo & Giacalone 2014). Here, we note the difference between our assumption on the 
waves and theirs. We assume propagating (superluminal) EM waves. On the other hand, 
their magnetic held is static in the huid rest frame (entropy waves), and the motional electric 
held E = —{v X B)/c is used to calculate the particle acceleration in the shock rest frame, 
where v and B are the huid velocity and magnetic held, respectively. This diherence has a 
big impact on the particle acceleration mechanisms such as the phase locking. 

The electric held and magnetic held are expressed by the superposition of elliptically- 
polarized waves. They are decomposed to linearly polarized sinusoidal waves. Strictly speak¬ 
ing, the wave form of SLSW does not have a pure sinusoidal shape, but rather a sawtooth 
like shape (Max & Parkins 1971). Our approximation of sinusoidal wave is adequate for 
ca 3> cUp. We assume that there are a primary wave which is generated from entropy mode 
and isotropically distributed daughter EM waves. It is assumed that the entropy mode is 
completely transformed to the primary wave for simplicity. The primary wave is assumed to 



be linearly polarized and propagates to 2 ;-direction as 


(9) 

( 10 ) 


Eo ^ /locos(tjol - 

Eo ^ (^o//9ph,o) cos(cjol - fco2)es,> 


where Aq, ojq, ko, /3ph,o are the amplitude, frequency, wavenumber and phase velocity, respec¬ 
tively. Here Sx and Sy are the unit vector toward a:-direction and ^/-direction, respectively. 
The secondary components are described as 


T^sec ( 2 ), t) 
k^sec k) 


N 


^H„exp{i(fc„ • f - Wnt + Cn)}kE,, 


n=l 

N 


A 


exp {i (4 ■X-Uj + (n)}kB,n, 

1 Pph,n 
n=l ^ ’ 


( 11 ) 

( 12 ) 


where An, kn, oOn, Cn are the amplitude, wavevector, frequency, and phase of each mode, 
respectively. Here /9ph,n = t’ph,n/c is the phase velocity of each mode, which is calculated 
from the dispersion relation (equation ([7])) as we will show later. Since the phase velocities 
are larger than unity, the amplitude of the electric held is larger than the magnetic held. 
We note that it can be understood by considering the Faraday’s law V x E = Each 

waves propagate toward = kn/\kn\, with which e'x,n and e'y^n form the orthogonal 
coordinate system. The polarization vector are written as ^E,n = cos'0ne'a;,n + * sin 
and ^B,n = —^ sin+ cos'0ne'y,n- Since the distribution of the secondary waves is 
uncertain, we assume it to be isotropic. In this paper, e(, „ is chosen randomly to make the 
distribution of the secondary waves isotropic. The polarization {^E,n and ^ 3 , 71 ) and phase 
are also randomly distributed. The amplitude of each mode is given by 

(13) 


c ^ a 

■^n ^sec^n 


N 


Eg. 

n=l 


where represents the mean intensity of the secondary waves. The number of Fourier 
components N is 10^ in this paper. We use the following form for the power spectrum 


1 + i^nTc)^ 


(14) 


where Tc is the coherence time which is set as oooTc = 1, and a — 2 is the power law index of 
the energy spectrum of the secondary waves. We set a = 11/3, which makes Kolmogorov-like 
turbulences. Here, AcUn is chosen such that there is an equal spacing in logarithmic cn-space, 
over the hnite interval cUmin A (x < cumax- We set the minimum frequency and the maximum 
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frequency as cUmin = <^0 and cumax = IO^Wq in all calculations. The sum of the EM energy 
density of the primary and secondary waves is set to be constant in each rnn as 




2 


_ 1 _ (-2 
^0 ‘'set 


(15) 


Using the nnit of time in this paper is dehned as 

me _ ^ 
eq 

The strength parameter is dehned by using ^ and loq as 

eq 

a = -. 

mojjQ 


(16) 


( 17 ) 


When there is no primary wave, the strength parameter is dehned by replacing ojq to cUmin- 
The phase velocities /5ph,n = {knc) are calculated from the dispersion relation 


u 


2 

n 


_2^ 

\/l + 


+ kte 


.2 „2 


(18) 


Here we note that the superposition of the modes is an approximation, because this dispersion 
relation is nonlinear. The appropriateness of this approximation for the obtained results 
should be checked in future works. We assume that the electrons initially have the relativistic 
energy yomc^ = lOmc^ and have isotropic velocity distribntion which mimics thermalized 
particles in the precursor region of the termination shock. In this paper we use 10^ electrons 
for the calculation. Since we neglect the back reaction from particles to EM helds, the 
electron nnmber is important only for the statistics and does not change the physics. The 
injection points of these particles are chosen to be homogeneous. We fix the plasma frequency 
as LJp = \/5 X lO^^cuo and it does not change in time. This fulhlls the stable propagation 
condition. It is lower than the plasma freqnency which is estimated from the initial thermal 
Lorentz factor. However, as the Lorentz factor of the electrons becomes higher, the plasma 
freqnency becomes lower. To avoid complexity, we approximate the plasma frequency by a 
hxed valne, which corresponds to 7th ~ 100. 

We inject relativistic electrons as test particles in this prescribed EM held. We solve 
the eqnation of motion 

^(ymeu) = -e(E + ^ x B) (19) 

by using the Buneman-Borris method. We neglect the radiation back reaction. This is an 
adequate approximation for the electrons which have not been injected in shock crossing 
process, because the cooling timescale is mnch longer than the dynamical timescale (Amano 
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& Kirk 2013). The radiation spectra are calculated from the information of the motion by 
direct using the Lienard-Wiechert potential 


dW 

dujdVL 


dvrc^ 


(1 — /5 • n)2 c 


( 20 ) 


where (3 = v/c\s the velocity of the electron, n is the observer direction, and t' is the retarded 
time (Hededal 2005). This formula can be applied for the frequency range of ca > yWp for 
the radiation in the plasma. In the next section, we will calculate the radiation spectra. We 
will see that the above condition satisfied almost all frequency range. When the exception 
is encountered, we will note it. 


3. Results 

3.1. Particle acceleration 

First, we demonstrate one of the notable features of the particle acceleration by the 
SLSW, namely the strong acceleration toward the wavevector direction. For this calculation, 
we set ujq = 0.1. Since the unit of frequency is e(;/mc as noticed earlier, the strength 
parameter is a = e<,/mcjjjQ = 10. We show the distribution of the x and 2 ; components of a 4- 
velocity at f = 3 x lO^cUf^^ in Fig[Tl The red cross’ are the 4-velocities for ec.sec/'^c = 1, which 
means that the turbulence consists of the secondary waves without primary wave. On the 
other hand, green dots are the 4-velocities for e<^sec/’^c = 0.1. The corresponding amplitude 
of the primary wave satisfies eEo/mc = — ^g^g ^,)/me ~ 0.995, and therefore the primary 

wave is dominant in this case. The red cross’ distribution is nearly isotropic, while the 
green dots’ distribution is quite anisotropic. The anisotropy is due to the acceleration of 
the electrons toward the wavevector of the primary wave. Even though there are other 
waves, the primary wave can dominate the electron motion and the phase locking occurs for 
e<^sec/’^c = 0.1. The distribution for x-direction is symmetric and the width is around 10 
times smaller than that for the case with eq^ec/^c = 1. This is because the amplitude of the 
secondary components <^sec is 10 times smaller. 

In Fig [21 we show the energy spectra for different amplitudes of secondary components 
{ec.sec/mc = 1, 0.5, 0.1, and 10“^) at t = 3 x lO'^Wf^^. For ew/’^c = 10“^, the energy 
distribution has a pure power law distribution. In this case, the secondary components are 
quite weak compared to the primary component. As a result, electrons which are initially 
moving nearly parallel to the wavevector of the primary wave are in the state of ’’phase 
locking” and selectively accelerated. The cutoff energy reaches an expected value of 0^70 = 
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10 ^ in the timescale of o^'JqUq^ = which is the typical oscillation timescale of 

a ’’phase locked” particle. The other electrons also tend to accelerate to the wavevector 
direction of the primary wave, but this acceleration is weaker, because they are not in the 
phase locking state. To clarify the selective acceleration, we show the time series of the 
energy spectra for ew/’^c = 10“^ in FigO The high energy cutoff evolves in time due to 
the selective acceleration. On the other hand, the peak Lorentz factor does not vary and 
remains at 7 ~ 10, since the acceleration by secondary waves is quite weak in this case. This 
distribution shows a power law shape of dN/dj oc 7 “^. We do not intend to claim that this 
power law index is universal, because it can be altered by initial conditions. For example, 
the initial velocity (direction) distribution, assumed to be isotropic, apparently changes the 
resultant energy distribution, since the phase locking effect is strongly dependent on the 
angle between the velocity and wavevector of the SLSW. Here we stress that we obtain high 
energy electrons which have an energy much higher than the peak energy, and the phase 
locking effect by the strong wave is an important factor to realize it. 

For ew/’^c = 1, many EM waves accelerate the electrons without long term phase 
locking, which is prevented by the disturbance of the velocity direction by the other EM 
waves. As a result, electrons diffuse in the momentum space, and the distribution function 
/(p) is well described by the Gaussian function. We can see that the energy distribution 
dN/d'y oc p‘^f{p) is consistent with the power law distribution in the low energy side with 
the index +2 and an exponential cutoff in the high energy side. 

Lastly let us focus on the cases for e^sec/'n^c = 0.1 and 0.5. In these cases, both the 
primary and secondary waves affect the particle distribution (green and light blue line in 
Fig [2]). While most particles are diffusively accelerated by the turbulence (secondary waves), 
a small fraction of particles are selectively accelerated by the primary wave. This selection 
can clearly be seen in 4-velocity distribution for eq^ec/f^c = 0.1 in FiglU The peak energy 
difference between two distributions comes from the energy density of the secondary waves. 
The electrons in this peak energy range are accelerated diffusively by the secondary waves. 
Therefore, the energy density of the secondary waves is smaller and the peak energy is lower. 
If the Ao is slightly larger than <^sec in the upstream of the termination shock, this acceleration 
mechanism can produce a broad energy distribution which contains a high energy power law 
tail. 


3.2. Radiation 

The radiation spectra of the electrons moving in the SLSWs for ew/'nrc = 0.1 are shown 
in FigslHandEi We calculate the radiation spectra in the restricted time from t = 3x lO^ca,)'^ 
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to t = (3 X 10^ + 2 X The starting time t = 3 x corresponds to the time of 

the energy distribution of Fig[2j This integration timescale is longer than Photon Formation 
Time of the typical frequency of the synchro-Compton radiation by a factor of 2 x 10^. This 
ensures that we can resolve the radiation spectrum down to ~ 10 ^ times lower frequencies 
than the peak one. On the other hand, the time step for pursuing the electron motion is 
10“^. This ensures that we can resolve the radiation spectrum up to ~ 10^ times higher 
frequencies than the peak one. Strictly speaking, the time step should be smaller than the 
inverse of the radiation frequency, and it is much shorter than 10“^. However, it is showed 
that the radiation spectrum is well described by using the large time step which is about only 
one order shorter than the photon formation time (Reville & Kirk 2010). The horizontal 
axis is the frequency normalized by eq/mc. The vertical axis is the flux in an arbitrary unit. 
The jaggy lines are the calculated radiation spectra, and the smooth lines are the analytical 
synchrotron curve, which are shown for comparison. 

In FiglU the observer direction is nearly along the z direction. To be precise, we set 
the observer direction n = {nx,ny,nz) = (10“^,0, \/l — 10“^). The reason for not choosing 
n = (0, 0,1) will be explained later. The peak frequency of the radiation spectrum in FigS] 
is cjpeak ~ 10^. This can be understood by using the synchro-Compton theory (Rees 1971, 
Gunn & Ostriker 1971). The cutoff Lorentz factor y^ut is around 10^ for e^sec/'^^c = 0.1 
at f = 3 X as is seen in Fig [2l The peak frequency can be estimated as cUpeak ~ 

7 ^ew/'^c = 10^, in the same manner as the synchrotron radiation. This can be justified as 
follows. First, we consider an electron which is in the phase locking state in a strong wave 
and tentatively ignore the other waves. In this case, the Lorentz force is negligible except 
for the direction parallel to the velocity. The resultant trajectory is nearly straight, and the 
curvature radius of the orbit is very long compared to ymc^/ew- Next we add the other 
waves which are isotropically distributed. The electron trajectory obtains a wiggling shape 
with the curvature radius for each case to be ~ 7 mc^/e<^sec- The deflection angle during a 
typical deflection is ~ I/ 7 , since the strength parameter defined by using only secondary 
components such as ew/'^cwmin is unity. As a result, electrons emit the ’’synchro-Compton” 
radiation around the peak frequency. The large part of this radiation power comes from the 
electrons moving nearly along the x-direction, since the selectively accelerated electrons are 
moving in this way. Thus, the radiation spectrum has a clear peak at w ~ ~ 10^. 

We note that the deviation from the exponential cutoff in the highest frequency region 
comes from the jitter radiation contribution, since there are secondary waves which have 
higher frequency than cjmin = e<^sec/'^c (cf. Teraki & Takahara 2011). The spectrum in the 
frequency region lower than the peak is harder than the “isotropic” synchrotron theoretical 
one. In our case, the velocity distribution is quite anisotropic. The spectral index of the 
synchrotron radiation toward some direction, due to an electron, is 2/3 (cf. Jackson 1999). 
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The calculated spectrum seems to be slightly harder than oc so that additional 

mechanisms may contribute to this spectrum, but we do not discuss this topic further. More 
detailed analyses will be done in our future work. The important thing which we should 
stress here is that the primary wave does not contribute to the radiation directly. It works 
only for the energy gain of the electrons. The radiation power and typical frequency are 
determined by the energy of the electrons and EM energy density of the secondary waves. 

In Fig El we set n = (1,0,0). The spectrum is well described by the isotropic synchro- 
Compton radiation. The spectral index at the frequency region lower than the peak coincides 
with 1/3, and the peak frequency cupeak ~ 10^ is the expected value. The Lorentz factor of the 
electrons moving around the x-direction is 7 ~ 30 (see Fig [1] and Fig [2]). The magnetic held 
(^/-direction) of the primary wave mainly contributes to the radiation power for this case. 
Thus, the peak frequency is described as 'j'^eBoy/mc ~ (30)^ x 1 ~ 10^. The spectrum at 
the frequency region higher than the peak shows slower decline than the exponential cutoff. 
This is due to the superposition of the contributions from electrons with higher energies. 
Jitter radiation components do not stand out in this spectrum, because of eBoy/mcoJo 3> 1 
(cf. Teraki & Takahara 2011). 

To conhrm the fact that the scatterers needed for the radiation of FiglHare the secondary 
components, we show the radiation spectra for ew/'^c = 10“^ in FigO The secondary 
waves are extremely weak compared to the primary wave, so that the primary wave causes 
the radiation. The integration time is identical to that used for depicting Figs 0] and El 
The lower (blue) line is the radiation spectrum for the observer located in the x-direction. 
The peak frequency at ~ 100 can be explained in the same manner as above. The By 
component of the primary wave B^y strongly deflects the electron, and produces the peak 
frequency at ~ 'y'^eBQy/mc ~ 100, since the typical Lorentz factor for the unselected electrons 
(which is not moving toward the ^-direction) is ~ 10 and eB^y/mc ~ 1. The spectral shape 
roughly coincides with the theoretical synchrotron curve as seen in Fig El On the other hand, 
different features are found in the spectrum for the observer located in fi = {nx,ny,nz) = 
(10“^, 0, y/l — 10“^), which is shown by the red line. It is noted that the spectrum is shifted 
vertically by a factor of 200 to see the shape clearly. The spectral shape is clearly different 
from the synchro-Compton radiation, because this radiation signature directly reflects the 
nonlinear orbit. This large-scale orbit is determined by the strong primary wave, and the 
small-scale deflection angle is much smaller than I/ 7 , since the strength parameter for the 
secondary waves ew/mca' m^n = 10“^ is very small. The sweeping behavior (of the beaming 
cone) is completely different from the sweeping with curvature radius 'jmc^/eB, which is 
assumed for the synchrotron radiation. The particle trajectory is nearly straight, and the 
sweeping time is much longer than mcjec, and I/wq- In this case, the sweeping timescale is 
roughly 10 times shorter than the maximum phase locked oscillation timescale 27ra^7ga;(/^. In 
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our calculation, the initial velocity directions do not coincide the wavevector of the primary 
wave, and \E\ ^ \B\. Moreover, the phase velocity of this wave is larger than c. Such effects 
shorten the sweeping timescale by a few times compared to the case for the vacuum. As 
a result, the typical frequency is 7 max'^o/( 7 o®^) ^ ^( 10 ) ~ 10 ^- We note that the second 
harmonics can be seen around uj ~ 1000. We note that this radiation spectra in the lowest 
frequency region a; ~ 10 is not precise because ycup ~ 10^ x \/5 x 10“^ is 0(1). However, 
this frequency range is not important for the current discussion. Here, we showed that 
the scatterer for the radiation for ew/'^c = 10“^ is the primary wave. Furthermore, we 
conhrmed that the radiation signature is not always described by the synchro-Compton 
theory. From this fact, we can understand that the scatterers which realize the spectra in 
Fig m are the secondary waves. 

Lastly let us explain the reason for choosing the observer direction to be f? = {rix, Uy^Uz) = 
( 10 “^, 0 , \/l — 10 “^), and not n = ( 0 , 0 , 1 ). As mentioned above, the selectively accelerated 
electrons tend to move toward the z-direction, but there are very few electrons moving ’’very 
close” to the z-axis. The reason is as follows. When the phase locking occurs, the perpendic¬ 
ular forces (to the velocity) from the electric magnetic helds nearly cancel each other. From 
the balance of perpendicular forces, we can estimate the angle 9 between the wavevector 
and the velocity. By using the assumed parameter = \/5 x 10“^, 6 ^ ~ I /7 is obtained for 
7 < 10^. On the other hand, the angle is constant as 6 ^ ~ 10“^ for 7 > 10^. The electrons 
moving along the z-axis are slightly deflected and tend to have angles > 10 “^ from the x- 
axis. Since we want to see the radiation from the strongly accelerated particles, we set the 
observer in the direction n = {rix, Uy, Uz) = ( 10 “^, 0 , \/l — 10 “"^). 


4. Discussion 

The superluminal strong waves (SLSWs) should exist around the termination shock of 
the pulsar wind nebulae. They may play an important role in the particle acceleration. Par¬ 
ticularly, the particle acceleration by the SLSWs may work as a pre-acceleration mechanism 
for the Diffusive Shock Acceleration (DSA). We note that for the electrons undergoing DSA, 
the contribution on the ’’energy change” is small since it is much slower than the Bohm 
limit of DSA. To make the injection rate higher, the energy distribution should be broader. 
The amplitude ratio between the primary and secondary waves is a key point for it. If the 
upstream EM held mainly consists of primary SLSW, the resultant electron energy distribu¬ 
tion tends to show a power law shape, and the 4-velocity distribution is anisotropic. On the 
other hand, if the SLSW distribution is isotropic and no primary wave exists, the obtained 
energy distribution is narrower than the former case and isotropic. The SLSW acceleration 
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in upstream may determine the injection rate of the DSA. The radiation from the electrons 
accelerated by the SLSW can be understood as the synchro-Compton radiation in the sec¬ 
ondary waves. In this section, we discuss the applicability of this acceleration to the PWNe 
and the observational prospect. 


4.1. Length scale of the acceleration region 

The length scale the SLSWs can exist is also unsolved, but is an important problem. 
Here we estimate the length scale needed for the acceleration. We consider the observer 
frame {K frame) and the upstream rest frame {K' frame). First, we consider the case for 
which these electrons are diffusively accelerated. We assume the velocity direction is nearly 
parallel to the boost direction of the K' frame in the K frame. The length scale in K' frame 
for a typical deflection is ~ c/a;Tjj, where is the typical frequency of the waves. This 
scale is written in the K frame as L ~ r(n' -|- H)/a;^in. In the K frame, we have 

L ~ IT^C/Wmin, (21) 

where Wmin = tuo is the frequency in the iF-frame which can be regarded as the inverse of 
the spin period of the pulsar. It is estimated by using the Crab parameters as 

i~6xlO‘TTf cm, (22) 

Next we consider the electrons in the phase locking state. The length scale in the K' frame 
is ~ a^ 7 oc/ci;o, where cJq is the frequency of the primary wave. In the K frame, it is 

L.6xl0-(^)'(|)hm. (23) 

This is shorter than the radius of the termination shock of Crab Nebula only by a factor of 
5. From this estimation, it can be realized that the maximum Lorentz factor reached in our 
simulation is an upper limit, since it may not oscillate even during in a period. On the other 
hand, the planer wave approximation becomes inadequate if the length scale is comparable 
to the termination shock radius. Therefore, the acceleration by secondary waves without 
phase locking can play a role in the upstream, but we should be careful about the phase 
locking acceleration. In other words, the maximum energy of the strong wave acceleration 
can be determined by the scale length the SLSWs exist. 
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4.2. Observational prospect 


Lastly we discuss the possibility of observing the signature of the pre-acceleration by 
the primary wave. The typical radiation in this situation is the synchro-Compton radiation 
from selectively accelerated electrons. The maximum frequency of the synchro-Compton 
radiation in the K' frame is written as 


^max 


r\j 


I max 


e<^ 

me 


(24) 


The maximum Lorentz factor may be determined by the scale limit, as we have discussed 
above. However, here we suppose the acceleration region is sufficiently large to reach the 
maximum energy in one cycle of the phase locking electron y^ax ~ ®^7o- We note that the 
Photon Formation Length of the typical synchro-Compton photon (cf. Teraki & Takahara 
2014) is much shorter than the length scale of the whole orbit. Therefore, we can estimate the 
typical radiation frequency without the information of the whole orbit. Here we suppose that 
the energy density of the SLSWs and electrons are highest near the region of the termination 
shock. The typical frequency of this radiation in the observer frame is estimated as 
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(25) 


Thus, it can be observed as a radio knot. The observable area is restricted as (rTs/r)2 by 
the beaming effect of the upstream bulk motion. Interestingly, small radio knots with scale 
< lO^^cm are observed near the termination shock by VLBI imaging (Lobanov, Horns, & 
Muxlow 2011). If we adjust the observable area to this observation, the bulk Lorentz factor 
is constrained to T > 10^. This value is consistent with our scenario and acceptable for 
the PWNe models. If these radio photons are emitted as in our model, unfortunately, no 
high energy counterparts emitted by the inverse Compton scattering will be observed. The 
luminosity should be much smaller than the gamma-ray from the whole nebula, and current 
gamma-ray observations cannot resolve spatial structure of the PWNe. 


5. Summary 

We have investigated the electron acceleration in the superluminal strong waves and the 
radiation from them. We considered two classes of waves. One is the primary wave, and the 
other is the isotropically distributed secondary waves. We took the amplitude ratio of them 
as a parameter. When the primary wave is dominant, the electrons moving nearly along the 
wavevector direction are selectively accelerated, and form a power law distribution in the 
high energy region. On the other hand, when the secondary waves are dominant, the energy 
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distribution is narrow and shows an exponential cutoff, and does not show any power law 
tails beyond the peak. We can expect both cases in the upstream of the termination shock of 
the PWNe. If the former case is realized, this acceleration mechanism may play a signihcant 
role for the injection to the shock acceleration. The radiation features can be described by 
synchro-Compton theory when the amplitude of the secondary waves are similar to that of 
the primary wave. The radiation from the phase locked electrons can show different spectrum 
when the primary wave is extremely dominant and the strength parameter of the secondary 
waves is much smaller than unity. However, such situation may not be realized around the 
termination shock of the PWNe, since it is a very coherent situation. Radio knots which 
are illuminated by synchro-Compton mechanism without counterparts in higher frequency 
range can be observed in the immediate upstream of the termination shocks if the SLSWs 
can propagate much longer than its wavelength. 
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Fig. 1.— Distribution of 4-velocity for e<^sec/’^c = 0.1 (green) and 1 (red) at f = 3 x lO^Wg 
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Fig. 2.— Particle energy distributions for t = 3 x and a = 11/3. The horizontal 

axis is Lorentz factor of the electrons, and the vertical axis is dN/d'y. The curved lines are 
corresponding to ew/’^c = 10“^ (red), 0.1 (green), 0.5 (light blue) and 1 (yellow). 
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Fig. 3.— The evolution of the energy spectrum for ew/’^c = 10“^. Each lines are the 
spectra for t = (red), (green) and 3000ci;|^^ (light blue). The green and light 

blue lines are shifted in vertical direction by a factor of 2 and 5, respectively, to see their 
difference clearly. 
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Fig. 4.— Radiation spectra from the accelerated electrons for = 0.1 at t = 3 x 

Vertical axis is the flux in arbitrary unit and horizontal axis is frequency in unit at 
eq/mc. The observer direction is fi = (10“^, 0, \/l — 10“^), which is nearly parallel to the 
wavevector of the primary wave e^. The curved light blue line is the synchrotron theoretical 
curve. The straight green line is F,^ oc and straight brown line is F,^ oc for 


comparison. 
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Fig. 5.— Radiation spectra from the accelerated electrons for ew/'^c = 0.1 at t = 
3 X The observer direction is in the x direction. The curved (green) line is the 

synchrotron theoretical curve for comparison. 
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Fig. 6. — Radiation spectra from the accelerated electrons for e<^sec/''^c = 10“^ at f = 
3 X The observer direction for npper (red)line is n = (10“^, 0, \/l — 10“^). The 

observer for the lower (bine) jaggy line is in the x direction. The cnrved (magenta) line near 
the npper line is the synchrotron theoretical curve. 





